Fourier's Law: insight from a simple derivation 
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The onset of Fourier's law in a one-dimensional quantum system is addressed via a simple model 
of weakly coupled quantum systems in contact with thermal baths at their edges. Using analytical 
arguments we show that the crossover from the ballistic (invalid Fourier's law) to diffusive (valid 
Fourier's law) regimes is characterized by a thermal length-scale, which is directly related to the 
profile of the local temperature. In the same vein, dephasing is shown to give rise to a classical 
Fourier's law, similarly to the onset of Ohm's law in mesoscopic conductors. 
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Introduction. - Fourier's law of heat conduction [1] 
states that when a system is subject to a temperature 
difference, a uniform temperature gradient Vr(r) ensues 
in its interior and the heat current density j (r) is propor- 
tional to that gradient, j(r) — ~kWT{r). The propor- 
tionality constant k is called thermal conductivity and 
could be position dependent. 

Although this law is almost two centuries old, a gen- 
eral demonstration, or a set of conditions for its validity 
in the general case, is still lacking [2, 3]. The problem 
has recently received renewed attention, both theoreti- 
cally [4, 5] and experimentally [6] , especially in quantum 
systems, due to the growing interest in energy transport 
at the nanoscale, whose understanding is an important 
step towards utilizing nanostructures for potential energy 
applications [7-9]. 

The difficulty in proving this law from first principles 
is twofold. For one, a local temperature (in an inher- 
ently out-of-equilibrium situation) needs to be defined. 
Secondly, a local heat current needs to be calculated and 
these two quantities need to be compared. Both of these 
quantities are difficult to evaluate for the simple reason 
that they are usually neither well defined (like the case of 
the local temperature [10]) nor their definition is unique 
(for the case of the heat current [11]). In the quantum 
case, an additional complication arises which is related 
to the fact that the size of the Hilbcrt space of a given 
system generally increases exponentially with system size 
[12], so that the problem easily becomes computationally 
very demanding. 

In recent years, several models have been put forward 
where Fourier's law has been demonstrated [13-26], all of 
which employ in some way the idea of local equilibrium. 
The local temperature is usually either calculated from 
the expectation value of some local energy operator [15, 
22, 27, 28], or a uniform temperature gradient is assumed 
to exist [11]. An alternative route is to study the energy 
diffusion in closed systems (i.e. without thermal baths) 
[29], or to study a system with self-consistent reservoirs 
[21, 24-26]. 

In Ref. [5] , the present authors utilized a novel scheme 
from the theory of open quantum systems [12] to cal- 



culate the local temperature and the onset of Fourier's 
law in electronic quantum wires. The main findings were 
that for a ballistic system Fourier's law is invalid, and the 
temperature is constant along the wire. A temperature 
gradient develops as disorder is introduced in the wire, 
eventually leading to the onset of Fourier's law. These 
results, however, relied on numerical simulations and no 
simple analytical form could be deduced. 

In this paper we aim at calculating analytically both 
the local temperature and the heat current for a model 
one-dimensional quantum system. Our results indicate 
that Fourier's law is related to a thermal length scale, ei- 
ther quantum (localization length) or classical (dephas- 
ing length), which is reflected in the local temperature 
profile. Fourier's law is then valid only if this thermal 
scale is smaller than the system length. In particular, 
dephasing gives rise to a classical Fourier's law, similarly 
to the onset of Ohm's law in mesoscopic conductors. 

Model. - We consider a set of N weakly-connected iden- 
tical sub-systems, S'„. Each sub-system 5'„ has its own 
Hamiltonian Hn = Efe efe|fc^"^)('^^"^l> where Jfc^")) are 
the many-body state vectors of S'„. The sub-systems 
are connected via a tunneling Hamiltonian ?i„.„+i = 

H.c. The full Hamiltonian is 

Here we assume that only the left-most and right-most 
sub-systems (S'l and S'r) are connected to external envi- 
ronments, and are consequently held at temperatures Tl 
and Tr, respectively. The other sub-systems are not cou- 
pled to an external environment, and hence transitions 
between states in the central sub-systems have to be me- 
diated by the edge sub-systems. The model is schemati- 
cally depicted in Fig. 1. 

Local temperature. - The aim of the calculation is to 
evaluate the local temperature of Sn as a function of 
Tl, Tr and n - which serves as a position variable. The 
main assumption of the model is that the tunneling in- 
teraction is weak enough such that each sub-system is in 
a local thermal equilibrium at a certain temperature T„, 
described below. 

Let us start by considering only T^l and a single ad- 
ditional sub-system Tii. By stating that S'l is main- 
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FIG. 1: (a) Schematic representation of the modeL The sys- 
tem consists of identicah weakly connected sub-systems. The 
left-and right-most sub-systems are held at temperatures Tl 
and Tr, respectively, (b) Representation of the transition 
process in Si (see text). 



tained at a temperature Tl, we assume that it is in 
a local thermal equilibrium. Due to the external envi- 
ronment (which determines the temperature), there are 
transitions between different states k and k' , which are 
defined via the scattering rates wj^^f,, . The statisti- 
cal meaning of the temperature is that there is a single 
number Tl, which characterizes all the different transi- 
tions by a single rule of detailed balance. This means 
that regardless of k and k' , the transition rates obey 

W^^^ , / Ac / \ 

y^'lV)''' = 6xp I rt*^ ) (taking ks = h = 1 throughout 

the paper), where Ackk' — Cfc' — ^k- 

Now consider the system Si. Since 5*1 is described in 
the energy basis, a transition between states in Si is an 
inelastic process, which (in the absence of interactions) 
requires the presence of an environment. Since the only 
subsystem in contact with an environment is Sl, for a 
particle in a state k in 5*1 to scatter to a different state 
k' in 5'i, a scattering event to some state fci in Sl has 
to first occur. From that state, another transition will 
occur to a state k2 in S'l and a final transition to k' (this 
process is depicted in Fig. 1(b)). Thus, the transition rate 
in 5i is given by 



.(1) 

• k^k' — ^ 

k\ ,k2 



(1) 



where F^^^, is a transition probability from the state 
in n to the state k' in n', and is therefore propor- 
tional to the overlap between jfe*^""^) and ^}, and 
hence to |Vfcfc'p (in similarity to the Fermi golden rule). 
We now assume that the main contribution comes from 
states of the same energy, and for simplicity take a uni- 
form tunneling Hamiltonian, i.e. Vkk — V, which yields 
kk' ■ This implies 



^ k~>k' — 



W,, 



(1) 



fc— fc' 



(2) 



where 7 = F^ (note that in our notation states la- 
beld by k may belong to different subsystems, denoted 



by the upper index (n)). Now, the temperature of Si 
may be obtained from the detailed balance of S'l, via 

'^'IT^''' = exp (— ^j"' ) . Employing Eq. (1), the 7 pref- 

actor cancels and we find that Ti = T^. This is a simple 
manifestation of the fact that two systems in contact with 
each other equilibrate. 

The next step is to consider a chain of N sub-systems, 
still only connected to a single Sl- The goal now is to 
find the temperature of Sn- One can repeat the procedure 
described above, with the only change being that n in- 
termediate, neighboring system transitions occur before 
the particles have a transition event at Sl- Therefore, 
a simple generalization gives Wl%, = 7"M^i^fe', which 
again gives T„ — T\^- This argument is valid to lowest or- 
der in 7, which physically corresponds to including only 
sequential tunneling processes. 

By considering the addition of the right-most system 
Ss, and the edge of the A^-long chain, one similarly ob- 
tains 



W, 



in) 



k^k' 



(3) 



In this case, due to the presence of the right temperature 
Tr one cannot simply cancel out the prefactor, and the 
expression for Tl becomes more complex. In order to 
make progress, we assume that Tr = Tl -I- (5T, ST « 
Tl, and assume that the transition rates take the form 
Wtl^ (X exp (-^^ 

tion that Sl and Sn are at equilibrium). Substituting 
back into Eq. (3) and taking the first order in 5T, we 
find (after some algebraic manipulation) 



cx exp ^— ^ (in agreement with the assump- 



T„ = Tr 



5T- 



1 



1+7- 



2n-N 



(4) 



Defining a position variable x = ^ — 1 we obtain 
^(^) = i+exp(a:/c) ' whcrc ^ = [N \og^)~'^ is the "thermal 
length" that defines the length-scale over which there is 
an (approximately) uniform temperature gradient. 

In Fig. 2 we plot (T(x) - Tl)/<5T (which scales from 
to 1) as a function of position variable x, for different 
values of^ — 0.1, 0.5, 1 and ^ = 00, which corresponds to 
7 = 1. For small values of 7 (and hence of ^), most of the 
temperature change is close to the middle of the wire. In 
those regions, a uniform temperature gradient is indeed 
developed, and hence Fourier's conjecture is valid. When 
7=1 the temperature is uniform in the wire, and hence 
no temperature gradient ensues. This figure should be 
compared with Fig. 1 of Ref. [5], which exhibit similar 
features, albeit obtained from a microscopic model (with 
no a priori assumptions). 

We note that while strictly speaking the case 7 = 1 is 
beyond the above perturbation analysis, this result is still 
valid. This is because in the 7 = 1 case all the wave func- 
tions are delocalized and span the entire system. Thus, 
S'l and Sr. have the same weight in the transition rates of 
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Sn regardless of n, giving rise to a uniform temperature 
(which is just Tn = Tl + ^ST for a small temperature 
difference). 
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FIG. 2: The (normalized) local temperature {T{x) - T-l)/ST 
as a function of the position variable x. This figure should be 
compared with Fig. 1 of Ref. [5] . 

Heat current. - The next step towards understanding 
Fourier's law in our model system is to calculate the heat 
current. If the sub-systems are weakly coupled to each 
other, the local energy is naturally defined [II] as En = 
(7i„). From the continuity equation [5, II] one has a 
simple expression for the heat current, j„ = —T{En+i — 
En) (taking the distance between the sub-systems to be 
a = I). This definition for the heat current could be 
understood by noting that (under the assumption of a 
constant T) TEn+i is simply the rate of energy flow from 
Sn+i to Sn, and vice versa, and thus j„ defined above 
describes the net energy flow between Sn+i and Sn per 
unit time, i.e., the heat current. 

To calculate j„, we assume that for each sub-system 
Sn we can define the probability P^"'' to find the system 

in the state j/c^"^). Then, we have En = X]/c ^'^'^ 
for the heat current 



(n+1) 



dn 



apj"^ dTn 
dTn dn 



(5) 



Since we are assuming local thermal equilibrium, it fol- 
lows that Pj.^^ oc exp(— efc/T„), and hence — — 

efc p(") 



This gives for the heat current 



Jn 



-rE 



in)' 



VT{n) ^ -K„Vr(n) 



(6) 



where k = ^2 i in agreement with the standard expecta- 
tions [30, 31]. Note that j is proportional to the tempera- 
ture gradient, and hence the thermal conductance scales 
as L~^, as required by Fourier's law. However, since 



the local temperature itself has a length-dependence, the 
above rule strictly applies only within the thermal length 
from the center of the wire. Since in experiments the 
measured k is a global property (that is, an average of 
K over the entire length of the sample) finite size effects 
may take place and give unusual scaling for k{L) [6]. 

The fact that this simple model exhibits Fourier's law 
at weak coupling can also be understood by comparing 
it to the results of Ref. [11], in which a general form of 
Fourier's law is derived for a general system which obeys 
three conditions. Our model satisfies these conditions, 
and hence one indeed expect it to display Fourier's law. 

However, we stress that the existence of a tempera- 
ture gradient does not necessarily imply "normal" heat 
conduction. While in a simple system as that described 
here we have shown that the two are connected, in more 
complicated systems, the non-linear nature of the inter- 
actions may give rise to a finite temperature gradient but 
anomalous conductance [32]. 

Dephasing. - Up to now we have assumed that there is 
contact between the chain and the external environment 
only at the edges of the chain. Let us discuss the effect 
of local environments acting along the chain, which may 
result in the dephasing of the wave- functions (in real sys- 
tems this may be caused by, e.g. inelastic scattering off 
low-energy phonons). We thus introduce a length scale 
L^, which characterizes the length over which the wave 
function retains its phase, li > N then the dephas- 
ing has no effect, and Fourier's law is valid in its form 
of Eq. (4). However, in the case < N, there is a 
new natural division of the system into subsystems of 
length L^. [33] Again, we assume that each of these sub- 
systems is in local equilibrium, and interacts weakly with 
the neighboring subsystem. However, as opposed to the 
case discussed above, there is no coherent tunneling be- 
tween the subsystems, but rather classical transport be- 
tween them. Thus, the rate equations for the occupation 
probabilities in the subsystems may be written as [34] 



where, for simplicity, we have assumed that the tran- 
sitions are between states of similar energies. Here 
n = I,2,...7V0 is the index of the sub-system, and the 
constant describes a typical transition rate between 
sub-systems defined via the division to sub-systems, 
and should in principal be determined microscopically. 
Equation (7), along with the boundary conditions 



P 



(L,R) 



- P 



(L.R) 



fc,cqilibrium 



OC exp(— ei;/TL,R,)), has a simple 



solution for its steady state, given by Pf 



P 



(L) 



^(R) 



J, ). Assuming that P^"' (x exp(— efe/r„) 



^^^^ 

and in the limit Tr — Ti^ = ST << Ti^ one obtains a lin- 
ear form for the local temperature, T„ « Tl + ;^<5P- We 
thus conclude that dephasing brings about the classical 
form of Fourier's law. Explicit examples can be found 
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in, e.g., Rcf. [21] and Rcf. [26], where local dcphasing 
was introduced (in a quantum harmonic lattice and elec- 
tronic system, respectively) by means of local external 
baths, giving rise to Fourier's law and a linear tempera- 
ture profile. 

It is also useful to consider the analogy between the 

effect of dephasing on the local temperature and on the 
resistance of a one-dimensional wire consisting of local- 
ized sub-systems (i.e. an Anderson insulator) [35]. In 
the absence of dephasing (L^ > N) the resistance is ex- 
ponential in the wire length, R ~ exp(A^/^ioc), where 
^loc is the localization length. However, when < N, 
the resistances of different subsystems of length (each 
with a resistance of i?^ ~ exp(L^/^ioc)) are connected in 
series, resulting in a linear dependence of the resistance 
on length, R ^ R^L/L^, i.e., classical Ohm's law. This 
crossover from a classical to a quantum regime is similar 
to what was demonstrated above in the case of Fourier's 
law. 

Discussion. - In summary, we have presented a sim- 
ple model where the local temperature and heat currents 
may be evaluated analytically. We have shown that the 
onset of Fourier's law requires the presence of a local ther- 
mal equilibrium at each sub-system that constitutes the 
full system. It breaks down in the case of strong coupling 
between the sub-systems. In that case, the temperature 
is constant throughout the sample. Including dephasing 
processes brings about a classical form of Fourier's law, 
in similarity to the onset of Ohm's law for the resistance. 

The generality of these results may be understood by 
the following arguments. Consider a one-dimensional 
system held at a temperature difference 5T. Now, the 
system may be broken into the smallest possible sub- 
systems {Sn in the above calculation). If each such sub- 
system has a unique temperature which defines the relax- 
ation rates between all the states in the sub-system, then 
the system is in local thermal equilibrium and Fourier's 
law is valid. If not, the system may be coarse-grained 
to generate larger sub-systems, and again for each sub- 
system a local temperature is sought. If local thermal 
equilibrium is eventually obtained, then the number of 
coarse-grained sub-systems describes the effective length 
of the system, and along with the interaction between 
the sub-systems, it gives the thermal length ^ which de- 
termines the length-scale for a uniform temperature gra- 
dient to ensue. However, if the coarse-graining procedure 
reaches the scale of the system, only a single temperature 
can be defined. In that case, the temperature is uniform 
across the sample, the thermal conductivity diverges and 
Fourier's law is invalid. A similar phenomenon occurs in 
the presence of dephasing, where the role of the (quan- 
tum) thermal length is played by the dephasing length. 

Recently, a crossover from ballistic to diffusive thermal 
transport was shown numerically to appear in a system 
with self-consistent reservoirs [25] . The crossover was de- 
termined from a length-scale related to the coupling be- 



tween the system and the baths, in a similar fashion to 
the dephasing length discussed above. Our results are 
in agreement with (and thus provide an intuitive explana- 
tion for) the numerical observations in Ref. 25. However, 
there are indications that in certain models there is no 
length-dependent crossover [32]. Therefore, whether the 
onset of Fourier's law is always a crossover phenomenon 
seems to depend on the specific model and a satisfactory 
answer to this question has yet to be found. 

This work has been supported by DOE under grant 
DE-FG02-05ER46204 and UC Labs. 
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